Stern- Gerlach Effect of Weak-Light Ultraslow Vector Solitons 
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We propose a scheme to exhibit Stern-Gerlach (SG) deflection of high-dimensional vector optical 
soliton (VOS) at weak-light level in a cold atomic gas via electromagnetically induced transparency. 
We show that the propagating velocity and generation power of such VOS can be reduced to 10 -6 c 
(c is light speed in vacuum) and lowered to magnitude of nanowatt, respectively. The stabilization 
of the VOS may be realized by using an optical lattice formed by a far-detuned laser field, and its 
trajectory can be deflected significantly by using a SG magnetic field. Deflection angle of the VOS 
can be of magnitude of 10 -3 rad when propagating several millimeters. Different from atomic SG 
deflection, deflection angle of the VOS can be distinct for different polarization components and can 
be manipulated in a controllable way. The results obtained can be described in terms of a SG effect 
for the VOS with quasispin and effective magnetic moment. 

PACS numbers: 42.65.Tg, 42.50.Gy 
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Stern-Gerlach (SG) effect, i.e., a particle with nonzero 
magnetic moment deflects when passing through an inho- 
mogeneous magnetic field, was firstly discovered in early 
time of quantum mechanics. This effect illustrates the 
necessity for a radical departure from classical mechan- 
ics, and characterizes quantum nature of atomic motion 
in a simple and fundamental way [lj. Recently, simi- 
lar effect was also predicted in many other systems, e.g., 
spinor Fermi and Bose gases [2J and chiral molecules [3|. 

All massive elementary particles, such as electrons, 
have non-zero magnetic moments. Contrarily, photons 
have no magnetic moment in vacuum, thus experience 
no force when passing through inhomogeneous magnetic 
field. Recently, in a very remarkable experiment 
Karpa and Weitz showed that photons may acquire effec- 
tive magnetic moments when propagating in a resonant 
atomic gas, and hence can deflect by a gradient magnetic 
field. In their experiment, a technique of electromag- 
netically induced transparency (EIT) |5| is exploited, by 
which a small absorption and slow propagating velocity 
of photons can be realized. 

However, the EIT-enhanced deflection of light in Ref. 
Q cannot be explained as a standard SG effect because 
only one component of "spin" is involved. In this Letter, 
we propose a double EIT scheme to demonstrate a SG 
effect of high-dimensional vector optical soliton (VOS), 
which not only has two polarization components (i.e. a 
quasispin) but also allows a distortionless propagation. 
Propagating velocity and generation power of the VOS 
can be reduced to very low level. Stabilization of the 
VOS can be realized using an optical lattice formed by a 
far-detuned laser field. The VOS can acquire very large 
effective magnetic moments, and the deflection of its tra- 
jectory is much more significant when passing though a 
SG gradient magnetic field. 

Before proceeding, we note that besides Ref. 3, opti- 
cal beam deflection in external fields has been the subject 
of many previous works The present work is re- 



lated to Ref s. 0, [13] and to recent studies of slow- light 
solitons [IlHU Essence of Refs. 0, S, EU is a SG ef- 
fect of linear polaritons. However, such linear polaritons 
spread and attenuate during propagation because of the 
existence of diffraction and other detrimental effects. In 
Refs. (ll| - [l3 |. slow-light solitons via EIT are suggested, 
but no SG effect is considered. 

In contrast, the scheme presented here exploits optical 
lattice and EIT-enhanced Kerr effect, which allow the for- 
mation and stable propagation of high-dimensional VOS, 
or called nonlinear polariton, with effective magnetic mo- 
ment (SG deflection) being four (two) orders of magni- 
tude larger than that of the linear polariton of Ref. 0] ■ 
Thus, comparing with that obtained in a linear scheme 
0> B E3] j the SG effect proposed here is more efficient 
and robust for observation and practical applications. 

To be specific, we consider a medium consisting of 
five-level atoms with M-configuration. A linearly polar- 
ized, pulsed probe field (with pulse duration to) E p 
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drives the transitions |3) «-» |2) and |3) f-> |4) by its left- 
circular (i.e. <r~) polarization component E p i and right- 
circular (i.e. er + ) polarization component E P 2, respec- 
tively. Here £ p i jP 2 are envelopes and e± = (x±iy)/V2- A 
7r-polarized, strong continuous-wave control field E c i = 
z£ c i exp[i(fc c ix — uj c it)) + c.c. (E C 2 = z£ C 2exp[i(k C 2X — 
0J C 2t)} + c.c.) drives the transition |1) O |2) (|5) «-» |4)) 
(Fig. [TJa) ). x, y, and z are unit vectors along coordinate 
axes x, y and z, respectively (Fig. (Ub)). 

We assume an inhomogeneous magnetic field B(y) = 
zB(y) = z(B + B x y) {B t < B ) is applied to the 
system. Here B contributes to a Zeeman level shift 
AEz = [isg'p'm'pBQ, and hence removes the degeneracy 
of ground-state sublevels \j) (j = 1, 3, 5) and the excited- 
state sublevels \l) (I — 2,4). ub, ff'p, and m F are Bohr 
magneton, gyromagnetic factor, and magnetic quantum 
number of the level \j), respectively. B\ contributes a 
transverse gradient of the magnetic field, resulting in a 
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FIG. 1: (color online) (a): Double EIT scheme. E p and E c j 
(j = 1, 2) arc probe and control fields, respectively; S p , 8 P + A, and 
5 c j are detunings. (b): Absorption spectrum Imffj (u) as functions 
of uj. Solid and dotted lines correspond to the a~ and cr+ po- 
larization components, respectively, (c): A possible experimental 
arrangement, where an inhomogeneous magnetic field B removes 
the degeneracy of ground states \j) (j = 1, 3, 5) and excited states 
\l) (I = 2,4), and causes Stern-Gerlach deflection of probe-field 
components. 8± and 82 are deflection angles of u~ polarization 
component (i.e. E p i) and t7 + polarization component (i.e. E p 2) of 
high-dimensional VOS, which has a quasispin and an effective mag- 
netic moment. The curved thick arrow represents the far-detuned 
optical lattice field E used to stabilize the VOS. 

SG deflection of the probe field. 

We assume further a small, far-detuned laser field 
E(x, t) — xE cos(x/R±) cos(wl^) is also applied into 
the medium, where Eq, R±, and u>l are field ampli- 
tude, beam radius, and angular frequency, respectively. 
Due to E(i,i), Stark level shift AE jjS = -\aj{E 2 ) t = 
— hctjE 2 (x) occurs, here aj is the scalar polarizability of 
the level \j), (•••)* denotes the time average in an oscil- 
lation cycle, and hence E{x) = (Eo/y/2) cos(x/R±). The 
aim of introducing the far-detuned laser field is to form an 
optical lattice potential to stabilize the high-dimensional 
VOS [Uj], as shown below. 

Besides, atoms are assumed prepared initially in the 
ground-state level |3) and trapped in a gas cell with ul- 
tracold temperature to cancel Doppler broadening and 
collisions. Thus, the system is composed of two A-type 
EIT configurations (i.e. double EIT). A possible arrange- 
ment of experimental apparatus is suggested in Fig. [He). 

Under electric-dipole and rotating-wave approxima- 
tions, the Hamiltonian of the system in interaction pic- 
ture is Htat/h = (S p - S el )\l)(l\ + 6 P \2)(2\ + (S p + 
A)|4)(4| + (5 P + A - <5 c2 )|5)(5| + fi c i|2)(l| + fi p i|2)(3| + 
r2j, 2 |4){3| +n c2 |4){5| +H.c, where Q p i=-(p 23 • e_)£ p i/fi 
and fip2=-(p 43 • e + )£ p2 /h (fl c i=-(p 21 • z)E c i/H and 
Q C 2=— (P45 • 2,)£ c2 /h) are respectively Rabi frequencies 
of two circularly polarized components of the probe field 
(two 7r-polarized control fields), with p j7 being the elec- 
tric dipole matrix element associated with the transi- 



tion from \j) to \l). The detunings are defined as S p — 
U23 + M23-B(y) - ±a 23 E(x) 2 - uj p , S cl = LO21 + V21 B(y) - 
\a 21 E(x) 2 -ui cU S c2 = oj 45 + p 45 B(y) - ±a 45 E(x) 2 -us c2 , 

and A = n± 2 B(y) - \a i2 E(x) 2 , where pbji = n B (ff , F m J F - 
g l F m F )/h, aiji = (aj-ai)/h, and uj jt = (Ej-Ei)/h, with 
Ej being the eigenenergy of the state \j). 

The motion of atoms is governed by the Bloch equation 
for density-matrix p, 

d P /dt = -i[H int ,p}/h-T(p), (1) 

where T(p) is relaxation matrix representing spontaneous 
emission and dephasing (see Supplementary Material). 
Electric-field evolution is controlled by Maxwell equation 

V 2 E - (l/c 2 )9 2 E/9t 2 = (l/e c 2 )d 2 P/dt 2 , (2) 

where P = A/"Tr(pp) is electric polarization with Af 
the atomic density. Under slowly varying envelope 
approximation, Eq. ([2]) reduces to [id z + (i/c)d t + 
cV 2 L /(2w p )]Opi iP 2 - k 32 , 34P23.43 = 0, where = 
d 2 /dx 2 + d 2 /dy 2 and k 32 , 34 = AT|p 32)34 • e T \ 2 uj p /(2he c) 
with eo the vacuum dielectric constant. 

Linear propagation of the probe field in the absence 
of diffraction can be obtained by taking fl p i, p2 as small 
quantities and £>i, Eo as zero. Then one has £l p j = 
Fj exp{i[Kj(ui)z — uit)]} (j = 1, 2) with A"i j2 (w) = lo/c + 
«3 2 ,3 4 (w — di^)/Di, 2 (linear dispersion relation). Here, 
Fj are constants, D\ 2 = |f2 c i iC2 | 2 — (ui — di^)(cu — d 2j4 ), 
di = (S p - S c i) - 1713/2, d 2 = S p - i(T 2 + 7 23 )/2, d 4 = 
(S p + A)-i(T 4 + 734 )/2, and d 5 = (S p + A - 5 c2 ) - i l35 /2 
with S p = cj 23 + p. 23 B - Lu p , A = p. i2 B 0l S c i = lo 21 + 
M21-B0 - w c i, and S c2 = w 45 + p i5 B - w c2 . and 7^ 
denote the spontaneous emission and dephasing rates of 
relevant states, respectively. 

The linear dispersion relation displays two branches. 
Fig. [TJb) shows the absorption spectrum of lmKj(uj) 
(j = 1,2) as a function of frequency uj. Parameters 
are chosen for a laser-cooled 85 Rb atomic gas with |1) = 
\5 2 S 1/2 ,F = 2,m F = -1), |2) = |5 2 P 1/2 ,F = 2,m F = 
-1), |3) = |5 2 S 1/2 ,F = l,m F = 0), |4) = |5 2 P 1/2 ,P = 
2,m F = 1), and |5) = |5 2 S 1/2 ,F = 2,m F = 1). De- 
cay rates are r 2 ~ r 4 ~ 6 MHz and 7 i 3 ~ 723 ~ 
73 4 ~ 7 35 ~ 50 Hz. Other parameters are taken as 
K32 ^ K34 = 1.0 x 10 9 cm _1 s _1 , n c i = il c2 = 1.0 x 10 7 
s _1 , S p = S c i = S c2 = 0, and Bq = 34.1 mG. The solid 
(dotted) line in the figure is for er~ (c + ) polarization 
component. We see that large and deep transparency 
windows in the absorption spectra of both polarization 
components (double EIT) appear. Using above parame- 
ters, group velocities of the both components (defined by 
V gj = Re(dK j /dtj)- 1 ) are given as 3.3 x 10^ 6 c. 

However, the linear solution is unstable due to the 
diffraction and other detrimental effects, which results 
in spreading and attenuation of the probe field during 
propagation, as demonstrated by Eq. (24) of Ref. [Toj | . 
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To solve this problem we use nonlinear effect to sup- 
press the spreading and attenuation. When including 
weak nonlinearity and diffraction, we obtain the follow- 
ing nonlinearly coupled, dimensionless equations, derived 
by using a standard method of multiple-scales (see the 
Supplementary Material): 
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+5l2,2l|u2,l| 2 )"l,2 + Vi,2(C) 57)^1,2 = 

where r=t/r , (£,r))=(x,y)/R x , 
and Uj (r, £, r ? )=[fi f , i /^( / 9 i )] e - iRe ^l' 
Lniff = u!pR±/c and [To are respectively the typ- 
ical diffraction length and Rabi frequency, Fj(pj) 
are normalized Gaussian functions (i.e. Fj = 
[l/(Po\/^)] 1/2 exp[-p|/(2p2)] with Pj = (z- V g3 t)/L uis 
and po a constant 10]), 311,12,21,22 = Wu, 12,21,22/1^22! 
are nonlinearity coefficients with Wu t 22 = 
-K32,34di, 5 (|di, 5 | 2 + |Oci,c2| 2 )/Pi,2|-Di, 2 | 2 ) and 

W 12 ,21 = -«32,34dl, 5 (M 5 ,l| 2 + \^c2,cl\ 2 )/(D h2 \D 2A \ 2 ) 

characterizing respectively self-phase and cross-phase 
modulations, and Aj=Im[if,-| w= o]i>Diff are small absorp- 
tion coefficients. 

Combined potentials in Eq. have the form 



V j (Z,r,)=M j r 1 +M jC os 2 (0, 



where .Mi, 2 = Mi t2 R±Bi and A/1,2 = N\ <2 Eq are con- 
tributions from the SG gradient magnetic field (pro- 
portional to B{) and the optical lattice field (propor 
respectively. 



tional to Eq 
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as Mi, 2= -^32,34(^1,5/^23 
and Wi, 2 = ^32, 34 (c?i, 5 a23,43 
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Mj and V,- are defined 

|^cl,c2| 2 Ail3,53)-^Diff/-Di,2 
^cl,c2| ai3,53)-^Diff/-Dl,2- 

When deriving Eq. ©, B\ and Eq are assumed as small 
quantities. Additionally, tq is also assumed to be large 
(e.g. t = 2.1 ps) so that the second-order dispersion 
(proportional to d 2 Uj/dr 2 ) is negligible. 

Fig- E](a) and (b) show results of numerical simulation 
for \m | 2 respectively at t = and < = 3r for a deep opti- 
cal lattice (Eq = 3.2 x f 4 V cm -1 ). The soliton obtained 
displays a single-peaked structure. The result for |it 2 | 2 is 
similar to |iti| 2 due to symmetry and hence not shown. 
The case for a shallower optical lattice (Eq = 2.3 x f 4 V 
cm -1 ) is also simulated, with the result plotted in panels 
(c) and (d) for t = and t = 3tq, respectively. We see 
that in this case a multiple-peaked soliton appears. In 
both simulations, 8 P = 1.0 x 10 6 s _1 , S c2 = 1.0 x 10 5 s _1 , 
and R± = 16 pm with other parameters the same with 
those in Fig. [TJ In addition, Uq = 6.8 x 10 6 s _1 , which al- 
lows enough nonlinearity to balance the diffraction. The 
typical diffraction length LDiff and nonlinearity length 
i Non i (= 1/(Uq \W22D are ~ 0.2 cm. Furthermore, Bi 
is chosen as zero, i.e., the SG gradient magnetic field is 
absent, thus no SG deflection occurs. 

The stability of the high-dimensional VOS is checked 
by adding a small random perturbation to the stationary 



FIG. 2: (color online) (a) and (b): Evolutions of |iti| 2 respec- 
tively at t = and t = 3rn for single-peaked VOS. (c) and (d): 
Evolutions of |«i| 2 respectively at t = and t = 3to for multiple- 
peaked VOS. SG gradient magnetic field is absent (i.e. B\ = 0). 
The stability of the VOS is achieved by the far-detuned optical 
lattice. Result for |«2| 2 is similar to «i| 2 thus not shown. 
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FIG. 3: (color online) SG effect of ultraslow VOS. (a) and (b): 
Symmetric deflection (on j/-axis) of |«i| 2 and |it2 1 2 when propagat- 
ing from z = 2Lrjiff to 2 = 8Lrjjff (corresponding respectively to 
the subfigure from left to right), respectively, (c): Asymmetric de- 
flection of |ii2 1 2 (|ui| 2 is the same as (a) thus not shown), (d), (e), 
(f): Corresponding evolution of linear polariton. (g): Deflection 
angles of the VOS as functions of z/I/rjiff f° r B\ = 0.7 mG/zim. 
The solid line with positive (negative) slope is the analytical result 
of 0i ($2) for the symmetric case. Dashed line is the analytical re- 
sult of 9'2 for the asymmetric case (9\ is the same as the symmetric 
case thus not shown) . Points labeled by "x" and "+" are center po- 
sitions of the VOS polarization components obtained numerically. 



solution obtained in imaginary time (Fig. [5] (a), (c)) and 
evolving the solution according to Eq. ([3]) in real time. 
We find that the soliton can indeed propagate stably for a 
long time (Fig. [2](b), (d)). We have also used a standard 
linear stability analysis (see the Supplementary Material) 
to confirm the stability of the high-dimensional VOS. 

Next we study VOS deflection by numerically simu- 
lating Eq. © with Bi ^ 0. Shown in Fig. [3] (a) 
and (b) are spatial distributions of |wi| 2 (panel (a)) and 
\u 2 \ 2 (panel (b)) in (x, y) -plane when the VOS propa- 



4 



gates from z = 2I/Diff to z = 8LDiff with group velocity 
V g i ~ V g2 = 3.2 x 10~ 6 c. In the simulation, B 1 = 0.7 
mG /zm _1 is chosen. We see that an obvious deflection 
of VOS trajectories occurs due to the existence of the SG 
gradient magnetic field. Additionally, two different po- 
larization components deflect symmetrically in +y and 
—y directions, similar to the SG deflection for atoms. 

The SG deflection of VOS components can be made 
asymmetric. To show this, we take ft C 2 = 0.9 x I0 7 s _1 
without changing other parameters, then (V g i,V g 2) — 
(3.2, 2.6) x 10~ 6 c. As a result, the trajectory of cr~ com- 
ponent keeps unchanged, whereas the trajectory of o~ + 
component changes as shown in Fig. |3]Jc). This is dif- 
ferent from atomic SG deflection, where trajectories are 
always symmetric for two different spin components. 

For comparison, in Fig. |3fd), (e), and (f) we present 
results of corresponding evolution for a linear polariton. 
One sees that the probe pulse spread rapidly. Thus the 
nonlinear effect is necessary for obtaining stable VOS and 
its robust SG deflection. 

Analytical VOS solutions of Eq. fijfy can be gained un- 
der some approximations: (i)The small absorption term 
—iAjUj is disregarded. (ii)Since in the presence of the 
SG gradient magnetic field the two polarization compo- 
nents of VOS separate each other after propagating some 
distance, the cross-phase-modulation terms can be ne- 
glected. (iii)The optical lattice is deep enough so that 
Vj can be approximated as M.jf] +JV}(1 — £ 2 ). Tak- 
ing Uj{T,t,rj) = WjiTyrfifaffiexpliAfjVgjT], where <^(£) 
is the normalized ground state of the eigenvalue prob- 
lem {d 2 /di 2 - 2A£C 2 )0j = 2E i (j) j with £ e = —y/P^/2, 
and integrating out the variable £, Eq. ([3]) reduces to 
[(i/v gj )d T + (l/2)d*\ Wj - M^/{2^-K Po )g jj \w^w j + 

(Aijij— y/Afj/2)wj = 0, which admits exact soliton solu- 
tions [la]. A single-soliton solution (see the Supplemen- 
tary Material) gives 

n pj = u a 3 [i/( Pa ^)] 1/2 {yfWil*?i^ 

xe -{s-v g] Tf/(2 P l) e -V^« 2 /V^ secha?i (5) 

where A 3 = (2 5 /^ / \p /\g jj I) 1 / 2 , <fj = M jVgj r(r, - 

M jV 2 gj r 2 /6), and 0, = {2fif j ) 1 /*(<n-M j l$ j T i /2) (j = 
1,2). We see that both VOS components are localized 
in three spatial and one temporal dimensions. Thus, 
(wi,ii2) can be considered as a vector light bullet. 

After passing the medium with length L, the center 
position of the jth polarization component of the VOS 
is at (x,y 3 ,z) — (0,MjL 2 R±/(2L 2 liB ),L), with the prop- 
agating velocity along the z- (y-) direction given by V g j 
(Vj = MjV 2 jR±t/T 2 ). As a result, the expected deflec- 
tion angle of the jth VOS component is 

3 = Vj/V gj = (L/VgAOirti/pyBx, (6) 

where r = R^/Ljag, p = hk p is photon momentum, 
fJ, S oij = Mj Vgj hkp is effective magnetic moment. With 



the data in Fig. |3l we obtain fj, soi 1)2 = ±7.6 x 10~ 20 J/T, 
which is four orders of magnitude larger than the effec- 
tive magnetic moment for linear polariton of Ref. 0). 
From Eq. ^ we see the deflection angle of the jth po- 
larization component of the VOS is proportional to the 
medium length L, the SG gradient magnetic field B±, and 
inversely proportional to the group velocity V g j . In a me- 
chanical viewpoint, the deflection of the jth component 
of the VOS is caused by the transverse magnetic force 
Fj = ji ao \jB\ and deflection angles can be expressed as 
8j = FjtintjT /pj with tintj = L/V g j being the interac- 
tion time. Due to untraslow propagating velocity, large 
deflection angles may be observed even for small L. 

Fig- (3Kg) shows deflection angles of the VOS as func- 
tions of z/Lum for Bi — 0.7 mG//im. The solid line of 
positive (negative) slope is the analytical result of 9\ (62) 
by Eq. © for the cr~ (a+) component with V g \ ~ V g 2 = 
3.3 x 10 _6 c (i.e. the symmetric case). Points labeled by 
"x" are center positions of the VOS components obtained 
numerically. We thus have (81,02) — (1.6,-1.6) x 10~ 3 
rad for z — 0.4 cm, which is two orders of magnitude 
larger than that for linear polariton obtained in Ref. 
|4j. The dashed line is the analytical result of 82 with 
(V g i,V g 2) = (3.2,2.6) x 10~ 6 c (i.e. the asymmetric case) 
and points labeled by "+" are numerical results (6% is 
the same as the symmetric case). In both cases analyti- 
cal results agree well with numerical ones. 

The generation power of the high-dimensional VOS 
predicted above can be estimated by using Poynting's 
vector [ijj], which is 3.5 nW calculated using the above 
parameters. Thus, very low input power is needed for 
generating the VOS in the present double EIT system. 

In conclusion, a scheme is proposed to exhibit SG de- 
flection of high-dimensional VOS via a double EIT. The 
VOS has ultraslow propagating velocity and extremely 
low generation power. The stabilization of the VOS can 
be realized by using an optical lattice, and its trajectory 
can be significantly deflected by a SG gradient magnetic 
field. The results obtained can be described in terms of 
a SG effect of the VOS with quasispin and effective mag- 
netic moments. We expect that such large and robust SG 
effect may have potential applications in magnetometery 
and quantum information processing. 
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